Scaling of impact fragmentation near the critical point 
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We investigated two-dimensional brittle fragmentation with a flat impact experimentally, focusing 
on the low impact energy region near the fragmentation-critical point. We found that the universality 
class of fragmentation transition disagreed with that of percolation. However, the weighted mean 
mass of the fragments could be scaled using the pseudo-control-parameter multiplicity. The data for 
highly fragmented samples included a cumulative fragment mass distribution that clearly obeyed a 
power-law. The exponent of this power-law was 0.5 and it was independent of sample size. The 
fragment mass distributions in this regime seemed to collapse into a unified scaling function using 
weighted mean fragment mass scaling. We also examined the behavior of higher order moments of 
the fragment mass distributions, and obtained multi-scaling exponents that agreed with those of 
the simple biased cascade model. 

PACS numbers: 46.50.+a, 62.20.Mk, 64.60.Ak 
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I. INTRODUCTION 



From asteroids to nuclei, fragmentation phenomena 
can be seen everywhere. Many scientists and engineers 
have been fascinated by the fragmentation process, and 
much effort has been made to understand fragmentation 
1]. In particular, impact fragmentation of brittle solids 
has been investigated by simulation and experimentally 
H> El S. 01 El- The results have shown that cumulative 
fragment mass (or size) distributions exhibit a power- 
law dependence, in which the exponent depends on the 
dimensionality of the fractured objects. These results 
do not depend on the details of fragmentation, such as 
the material the fractured object is made from or the 
fragmentation method. Consequently, Oddershede et al. 
concluded that brittle fragmentation is a self-organized 
critical phenomenon 0, Q . 

Cumulative distributions obey power-law dependence 
when the imparted energy is sufficiently large, but brit- 
tle solids do not break when the imparted energy is very 
small. Therefore, one might assume that there is a tran- 
sition in the fragmentation phenomenon. Recently, Kun 
and Herrmann examined critical behavior in fragmenta- 
tion transition and concluded that fragmentation tran- 
sition belongs to the percolation transition universality 
class |8|. In addition, Campi has shown that nuclei frag- 
mentation has the same statistical properties as perco- 
lation By contrast, Astrom et al. performed other 
fragmentation model simulations and obtained a scaling 
law whose critical exponents differ from those of perco- 
lation |lCj . 

While there are some simulation results, there have 
been no experiments on critical fragmentation. For a 
proper analysis, the existing simulation results must be 
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compared to experimental data. Here, we experimentally 
investigate brittle fragmentation as a kind of critical phe- 
nomenon, and we report experimental results for glass 
tube fragmentation near the critical point. 



II. EXPERIMENTAL PROCEDURE 

Figure^is a schematic illustration of the experimental 
apparatus. We used glass tubes (50 mm outside diame- 
ter, 2 mm thick, and 150 mm, 100 mm, or 50 mm long) 
as the fractured objects. They have 2-D geometry, and 
it is easy to impart the impact vertically because they 
can stand by themselves. Kun and Herrmann Q and 
Astrom et al. [l(J | performed numerical simulation of 2-D 
samples. 

A glass tube was put in a plastic bag, which was in 
turn placed between a hard stainless steel stage (dimen- 
sions 400 mm x 500 mm x 50 mm) and a stainless steel 
plate (10 mm thick). The stage was fixed by placing 
heavy weights on it. A cylindrical brass weight with a 
flat bottom (3.77 kg) was dropped on the stainless steel 
plate vertically, and a planar failure wave propagated to 
the glass tube from the impact circle (cross section of the 
glass tube). Consequently, the glass tube was cleaved by 
the impact. This is a type of 2-D fragmentation. Four 
poles (10 mm diameter) were used as guide poles for the 
falling brass weight. Stoppers for the stainless steel plate 
were placed at each guide pole to prevent a secondary 
impact. While the samples were not annealed after the 
cutting processing, they had almost the flat end. We 
used a level in order to check the verticality between the 
weight and the stainless steel plate at each fragmentation. 
The variable in this experiment was the height that the 
weight was dropped from. 

After fragmentation, we collected the fragments in the 
plastic bag and measured the mass of each fragment with 
an electronic balance. The data were analyzed, and cu- 
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FIG. 1: Schematic illustration of the experimental apparatus. 
The weight and guide poles were made of brass. Other parts 
were made of stainless steel or steel, h is the falling height of 
the weight. 



mulative fragment mass distributions were obtained. We 
fractured a total of 60 samples (34 : 150 mm, 15 : 100 
mm, and 11 : 50 mm length). We defined the falling 
height ft, as a length from the top of stainless steel plate 
to the bottom of the falling weight. In our experiments, 
it was set in the range 61 mm < h < 407 mm. This 
range corresponds to 0.050 < e < 0.251 (Nm/g) in terms 
of the imparted energy (released potential energy of the 
dropped weight) per unit sample mass, e. 



III. RESULTS 

When the height was too low to break the glass tube, 
no visible cracks were seen. Once a visible crack was pro- 
duced, fragmentation occurred suddenly. We defined the 
point at which a sample began to cleave as the fragmen- 
tation critical point. We could not obtain samples that 
had only visible macro-cracks, but did not fragment. We 
show some photos of fragments in Fig. [21 We observed 
that with a flat impact, low impact energy fragmentation 
events produced vertical main cracks. Therefore, vertical 
cleaving produces a few, large fragments (Fig. |2{a)). In 
the mediate range, some large fragments remained (Fig. 
|2tb)). For large impact energy events, the largest frag- 
ment was smaller, and it was impossible to distinguish 
the main cracks from the many small fragments by ob- 
servation (Fig. IHc)). 








FIG. 2: Typical fragments photos resulting from: (a) small 
imparted energy, (b) intermediate region, and (c) large im- 
parted energy. Imparted energies per unit mass e for (a), (b), 
and (c) are 0.098, 0.118, and 0.139 (Nm/g), respectively. 



A. Cumulative distribution 

First, we plotted the cumulative mass distribution of 
each fragmentation. Figure|3]shows some typical cumula- 
tive fragment mass distributions, iV(m), for 150-mm-long 



samples. N(m) is defined as follows: 

/>oo 

N(rri) = / n(m')dm' , (1) 

J m 

where m and n(m) are the fragment mass and the num- 
ber of fragments of mass m, respectively. The distribu- 
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FIG. 3: Some typical cumulative fragment mass distributions 
for 150-mm-long samples. Power-law behavior is confirmed 
in the sufficiently high impact energy curves, while the dis- 
tribution curves for low impact energy do not show a clear 
trend. The corresponding values of fi for plus, diamond, tri- 
angle, square, and circle marks are 7.05 x 10" 4 , 3.32 x 1(T 3 , 
6.47 x ICT 3 , 2.34 x 1(T 2 , and 2.54 x 1CT 2 , respectively. 



tion curves for samples of different sizes were similar to 
the example in Fig. [3] One can see the clear power-law 
dependence for fully fragmented samples (far from the 
fragmentation transition point). The distribution curves 
near the critical point do not show clear power-law depen- 
dence, and they have a rather flat part. This tendency is 
quite different from that of percolation universality. In 
percolation universality, the clearest power-law depen- 
dence should appear beneath the critical point. In the 
fully fragmented data curves, we can observe a region 
that satisfies N(m) ~ to~( t ~ 1 ', i.e., n(m) ~ mT T . The 
obtained value of the characteristic exponent r is nearly 
1.5 (< 2). This result conflicts with the scaling ansatz 
of percolation Therefore, we believe fragmentation 
criticality to differ from that of percolation. Conversely, 
our result r — 1 ~ 0.5 concurs with Astrom's result |l0j . 
In addition, the value of r obtained in this experiment is 
consistent with Hayakawa's scaling for 2-D fragmentation 
with a planar failure wave Q- 



Since we are interested in fragmentation criticality, the 
height of the falling weight was controlled near the crit- 
ical point. Despite our control of the height (imparted 
energy), the M& data fluctuate too much due to individ- 
ual differences in the glass tube samples, e.g., the initial 
density of micro-cracks and residual stress distributions. 
When we considered the imparted energy as a control pa- 
rameter, we could not find clear critical behavior quanti- 
tatively. For example, we show a log-log plot of jj^ vs. e 
in Fig. 0fa). Where, is the weighted mean fragment 
mass. It is hard to say the data in Fig. 0fa) are par- 
ticularly convincing about showing linear behavior. We 
could not consider that power-law fitting was appropri- 
ate with confidence. Alternatively, we used multiplicity 
/i analysis, which Campi introduced to the analysis of 
nuclei fragmentation |9j- Multiplicity is defined as 



Mo 



(3) 



where Mo, Mi, and m m i n correspond to the number 
of fragments, the total mass of the fragments, and the 
smallest cut-off fragment mass, respectively. We fix 
f«min = O.Olg for all the analyses. This value usually 
corresponds to the smallest 2-dimensionality for the sam- 
ples. According to this definition, the fragmentation crit- 
ical point corresponds to fi = 0. While the multiplicity is 
a resultant parameter, rather than a control parameter, 
it can be considered the control parameter that incorpo- 
rates individual differences in the samples. Furthermore, 
because the multiplicity is easy to measure and calculate, 
it has been applied to other fragmentation systems 0. 
We use \x as a pseudo-control-paramctcr in this paper. 
Certainly, \i and the imparted energy are positively cor- 
related, roughly. The larger the impact energy imparted, 
the larger fj, becomes. 

In Fig. 01b), we plot |^ vs. of all the fragmenta- 
tion event data. It indicates the tendency toward di- 
vergence near the critical point. The criticality appears 
independent of sample size. All the data match the same 
power-law line. Hence, we assume that the relationship 
between the weighted mean fragment mass ^ and the 
multiplicity, fi, holds as, 



M2 

Mi 



(4) 



B. The weighted mean mass scaling 



The obtained value of a is 0.84 ± 0.05. 
nontrivial because a 1. 



This value is 



Since fragmentation events near the fragmentation 
critical point have less fragments, they cannot be ana- 
lyzed statistically using only N(m). Therefore, we used 
the fc-th order moment of the fragment mass distribution, 
Mk, defined as 

M fe = ^2m k n(m). (2) 



C. Scaling function 

We discuss the possibility of scaling N(m) curves using 
a. For large n data, N(m) curves have the same power- 
law, in which the power is r — 1 ~ 0.5. This value is 
independent of the size of the samples. This similarity 
allows us to collapse the N(m) curves of different-sized 
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FIG. 4: Critical behaviour of the weighted mean fragment 
mass ||^: (a) j^- vs. e. (b) vs. fi. The solid line indicates 
the form of the power-law Ma „ ( CT = 0.84 ± 0.05). 

samples into a unified scaling function. We expect the 
following scaling for distribution curves: 

P{m) ~ f(m/fi- a ) ~ m- (T ~ 1) g(m/^-' J ). (5) 

Where P(m) is the probability distribution of fragments 
whose masses are larger than m. We show the scaling 
result in Fig. [3] with the approximation P(m) = N ^™^ ■ 
The curves in Fig. include the results for samples of 
three different sizes (11 curves of large /j, data). Figure 
[3] shows good collapse of N(m) curves into the scaling 
function / and g. The function / has a clear power-law 
regime, in which the power is 0.5, and a unified cut-off 
scale. The former corresponds to the scaling region and 
the latter corresponds to exponential decay due to the 
finite size effect. Accordingly, the scaling function g has a 
constant part and a rapidly decaying part. The flat part 
in g exists only for the fully fragmented data, because 



FIG. 5: Plots of the scaling functions / and g approxi- 
mated using (a) N(m)/M ~ /(m/^~ CT ) and (b) N(m)/M ~ 
m ~( T ~ 1 '> g(rnj /i~ CT ). The values of r and a are obtained from 
scaling Figs. |H| and pjfb) as 1.5 and 0.84, respectively. 



clear power-law form N(m) ~ vnT^" 1 ^ can be satisfied 
with well fractured samples. The form of scaling function 
might change with /x, however, we consider they will not 
depend on the sample size. We could not confirm the 
clear collapse for small /i. More detailed experiments are 
necessary to verify the scaling functions in wide /i range. 

Bak et al. reported a similar scaling function for earth- 
quake statistics [12|. While they discussed the time cor- 
relation of earthquakes, their scaling function is analo- 
gous to ours. They concluded that the unified scaling 
function for earthquakes depends only on the number of 
earthquakes occurring within the area and period con- 
sidered. In Eq. |S1 our scaling function depends on the 
variable x — which is the normalized mass of the 

fragments. Moreover, the scaling Eq. |3]resembles the 
scaling function Astrom et al. obtained jlOj . 
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nature. However, this only alters the problem. The origin 
of the symmetry breaking given with a and the relation 
between a (or 07.) and r remain unsolved. The cascade 
model does not seem to yield the appropriate value of 
t (~ 1.5) directly. Since this model is too simple, it 
needs to be improved using more detailed analyses and 
experiments. 

Similar multiplicative model was proposed for the en- 
ergy cascade of eddy by Meneveau and Sreenivasan 0] . 
According to their result, the cascade of division into 
0.3 and 0.7 can explain the energy cascade of turbulent 
flow. However, the origin of the value 0.3 has remained 
unsolved. The value 0.3 slightly agrees the value -| we 
obtained. 



IV. DISCUSSION 



FIG. 6: Higher-order weighted mean fragment mass scaling. 
The marks indicate the experimental results, and the solid 
curve indicates the result calculated using Eq. Q(a = §)■ 



D. Multi-scaling and simple biased cascade model 

Next, we discuss the behavior of the higher-order mo- 
ments of the fragment mass distributions. We can easily 
expand Eq. ^ and introduce the higher-order weighted 
mean fragment mass scaling as, 



M, 



k+l 



(6) 



The open circles in Fig. show the Ofc values obtained 
from fitting the all experimental data like Fig. E£b). 
While the value (Jk=o is the trivial value 1, <Jk->oc ap- 
proaches a nontrivial value of approximately 0.6. 

In order to understand this behavior, we consider a 
simple biased cascade model. Let us set a unit mass 
fragment initially. First, this unit mass is divided into 
two pieces of mass a and 1 — a (0 < a < 1). The same 
a-biased partition occurs for each fragment in the next 
step. This cascade continues until the imparted energy 
dissipates. From the above definition and Eq. El the fol- 
lowing equation holds at each step: 



,fc+i 



+ (l-o) 



fc+i 



(7) 



Substituting a — \ into Eq. Q the trivial value o~k = 1 
is obtained for all fc. This case corresponds to mono- 
scaling critical fragmentation. In Fig. however, 
varies with k. This indicates a multi-scaling property of 
critical impact fragmentation. The solid curve in Fig. 
is the value calculated using Eq. [7| with a = 3 . The 
curve and experimental data roughly agree within the 
error bars. Although this model is very simple, it explains 
the trend of the higher-order behavior. We conclude that 
fragmentation near the critical point has a multi-scaling 



Kun and Herrmann have analyzed the largest fragment 
mass as an order parameter [SJ. In our experiments, the 
largest fragment mass had a very wide distribution, as 
seen in Fig. [5] while the weighted mean fragment mass 
behaved more calmly. Figure |21 suggests that the power- 
law exponent might change with fi. Ching et al. re- 
ported the impact energy dependence of the power-law 
exponent |l4j . Their interest focused on the very large 
imparted energy region and not on the neighborhood of 
the fragmentation critical point. In addition, if we fit the 
power-law form to all N(m) curves, the results involve 
large uncertainty, and do not show a clear trend. Ishii 
and Matsuhita reported that the distribution function for 
the low impact energy region has a log- normal form 0. 
We have not investigated the distribution of the concrete 
function form. The relationship between the two energy 
regions and the study of the concrete function form re- 
main open to study. 

Our experiments were focused on the 2-D fragmenta- 
tion with the flat impact. Kun and Herrmann inves- 
tigated the point impact 2-D fragmentation |8|. This 
might be a reason of the discrepancy between our re- 
sults and their ones. Hayakawa suggested that the uni- 
versality of the fragmentation depends not only on the 
dimensionality of fractured object, but also on the prop- 
agation dimensionality of the failure wave Q. Astrom 
et al. performed numerical simulations of 2-D fragmenta- 
tion with the flat expansion [T^j ■ Their results resemble 
ours in many points as described above. We believe that 
our results are valid for 2-D fragmentation by the flat 
impact to one side. Other dimensional objects or other 
failure wave propagation manner may produce other val- 
ues of scaling exponents such as r and o~k- Moreover, 
scaling-law among these exponents may exist. As writ- 
ten in Sec. [IJ the exponent r is independent of the ma- 
terial of fractured objects in fully fragmented state. We 
guess the other critical exponents are also independent 
of that, however, it should be examined. These problems 
are open questions. 

In general, the low impact energy region is difficult to 
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study. Small differences of initial condition are enlarged 
by nonlinearity of fragmentation dynamics. Figure E{a) 
presents this difficulty well. Even in Fig. Efb), the fluc- 
tuating data can be seen. Some unrealized experimental 
parameters might not be controled. We think that the co- 
incidence between the experimental result and the model 
shown in Fig. [S] is rough due to this reason. 

In conclusion, we investigated the 2-D impact frag- 
mentation of brittle solids experimentally. The measured 
characteristic exponents imply that the impact fragmen- 
tation of brittle solids does not belong to the percolation 
transition universality class. Instead, we found that the 
weighted mean fragment mass and the multiplicity 
parameter /i are related to the nontrivial scaling expo- 
nent a. The cumulative fragment mass distributions of 



different-sized samples could be collapsed into the uni- 
fied scaling function in the large \x regime. We calcu- 
lated the generalized weighted mean fragment mass scal- 
ing M J^ 1 ~ (i~ (7k . The behavior of the multi-scaling 
exponent was modeled using a simple biased cascade 
partition model. The critical behavior of other fragmen- 
tation systems should be studied in order to understand 
the universality of fragmentation phenomena in detail. 
Experiments with other materials, such as ceramics, or 
with samples of different dimensions may prove interest- 
ing. 

We thank Professor H. Sakaguchi and Professor K. No- 
mura for their useful comments. 
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